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Abstract: This paper summarizes some highlights from the Pierre Auger Observatory that were presented at the ICRC 
2011 in Beijing. The cumulative exposure has grown by more than 60% since the previous ICRC to above 21 000 
km 2 sryr. Besides giving important updates on the energy spectrum, mass composition, arrival directions, and photon- 
and neutrino upper limits, we present first measurements of the energy spectrum down to 3 ■ 10 17 eV, first distributions of 
the shower maximum, Jf max , together with new surface detector related observables sensitive to Jf max , and we present 
first measurements of the p-air cross section at ~ 10 eV. Serendipity observations such as of atmospheric phenomena 
showing time evolutions of elves extend the breadth of the astrophysics research program. 

Keywords: Pierre Auger Observatory, UHECR, cosmic rays 
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1 Introduction and Status of the Observa- 
tory 

The Pierre Auger Observatory started collecting data in 
2004. Since the completion of the base-line Observatory 
in 2008 its aperture has grown by about 7000 km 2 sr, each 
year. At this meeting we present data based on an exposure 
of more than 2 1 000 km 2 sr yr. The Auger Observatory uses 
hybrid measurements of air showers recorded by an array 
of 1660 water Cherenkov surface stations covering an area 
of 3000km 2 , together with 27 air fluorescence telescopes 
that observe the development of air showers in the atmo- 
sphere above the array during dark nights. 

An infill array with half the grid size has been completed 
and we present first data at this meeting extending the en- 
ergy spectrum down to 3 • 10 17 eV, thereby covering the 
ankle of the primary energy spectrum with full detection 
efficiency. Moreover, the three high-elevation telescopes 
(HEAT) started operation and - together with the infill ar- 
ray in the FOV of the telescopes - will allow us to extend 
the hybrid measurements further down to 10 17 eV with un- 
precedented precision. This will enable the study of the 
transition from galactic to extra-galactic cosmic rays. Con- 
struction of the buried muon detectors (AMIGA) in the in- 
fill area is in progress. Measurements of the muons are im- 
portant for studying the composition of cosmic rays from 
surface detector data and their information is also of vi- 
tal importance for studying particle interactions down to 
the energy of the LHC. In addition, an extensive R&D 
program for radio detection of UHE air showers is under 



way and construction of the Auger Engineering Radio Ar- 
ray (AERA) has started. Once completed, it will com- 
prise 160 radio antennas distributed over an area of 20 km 2 . 
First triple hybrid events composed of particle densities 
at ground, longitudinal shower profiles from fluorescence 
telescopes, and radio signals from the first antennas have 
already been observed. Last but not least, an intense R&D 
program for microwave detection of air showers has be- 
gun with the first GHz-antennas presently being installed. 
These extensions and new technologies may enhance the 
performance and capabilities of the Auger Observatory in 
Argentina and, in parallel, will explore their potential for a 
future much larger ground based observatory. 

There are 38 papers presented on behalf of the Pierre Auger 
Collaboration at this meeting, and they are all accessible in 
five e-print compilations with arXiv numbers 1107.4804, 
1107.4805, 1107.4806, 1107.4807, and 1107.4809. 



2 The energy spectrum 

An accurate measurement of the cosmic ray flux above 
10 17 eV is crucial for discriminating between different 
models describing the transition between galactic and ex- 
tragalactic cosmic rays, the suppression induced by the cos- 
mic ray propagation, and the features of the injection spec- 
trum at the sources. Two complementary techniques are 
employed at the Pierre Auger Observatory: a surface de- 
tector array (SD) and a fluorescence detector (FD). The 
energy spectrum at energies greater than 3 ■ 10 18 eV has 
been derived using data from the SD-array. The analysis 
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Figure 1 : Combined energy spectrum fitted with two power laws 
in the ankle region and a smoothly changing function at higher 
energies. Only statistical uncertainties are shown. The systematic 
uncertainty in the energy scale is 22 %. 



of air showers measured with the FD which also triggered 
at least one station of the surface detector array (i.e. hybrid 
events) enables measurements to be extended to lower en- 
ergies. Despite the limited number of events, due to the 
fluorescence detector on-time, the lower energy threshold 
and the good energy resolution of hybrid events allow us 
to measure the flux of cosmic rays with the standard array 
down to 10 18 eV, into the energy region where the transition 
between galactic and extragalactic cosmic rays is expected. 

The energy calibration of the SD-array is based on so- 
called golden hybrid events, i.e. events that can be indepen- 
dently reconstructed from the SD and FD data. Applying 
high quality cuts, 839 events could be used for the SD cal- 
ibration [1, R. Pesce]. The overall FD energy resolution is 
7.6 % and it is almost constant with energy. The total sys- 
tematic uncertainty on the FD energy scale is about 22%. It 
includes contributions from the absolute fluorescence yield 
(14%), calibration of the fluorescence telescopes (9.5 %), 
the invisible energy correction (4 %), systematics in the re- 
construction method used to calculate the shower longitu- 
dinal profile (10%), and atmospheric effects (6% - 8 %). 
The atmospheric uncertainties include those related to the 
measurements of aerosol optical depth (5 % - 7.5 %), phase 
function (1 %) and wavelength dependence (0.5 %), the at- 
mosphere variability (1 %) and the residual uncertainties on 
the estimation of pressure, temperature and humidity de- 
pendence of the fluorescence yield (1.5 %). 

The energy spectrum derived from hybrid data has been 
combined with the one obtained from surface detector data 
using a maximum likelihood method and is shown in Fig. 1 
together with a broken power law and a smooth cut-off at 
higher energies [1, F. Salamida]. Both, the ankle and sup- 
pression of the flux at higher energies are clearly visible. 
The spectrum can be compared to astrophysical models and 
can be described by both a proton and heavy-dominated 
composition at the highest energies. Thus, measurements 
of the composition are needed to discriminate between var- 
ious astrophysical models. 
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Figure 2: Centered distribution, X max — (X max ), for the lowest 
and highest energy bins. Subtraction of the mean allows the com- 
parison of the shapes of these distributions with the superimposed 
MC simulations. 



Data of the 750 m infill array reach full efficiency for all 
primaries at E > 3 ■ 10 17 eV and, using data with an expo- 
sure of 26 km 2 sr yr, extend the spectrum of Fig. 1 smoothly 
down to this threshold [1,1. MarisJ. Analysis of the com- 
position in this energy range is on-going. HEAT data, com- 
bined with the infill array, extend the energy range further 
down to 10 17 eV [2, H.J. Mathes]. 

3 The cosmic ray mass composition 

As discussed above, measuring mass composition of cos- 
mic rays along with the flux is a key to separating the differ- 
ent scenarios of origin and propagation of cosmic rays. The 
composition must be inferred from measurements of vari- 
ous shower observables, most importantly the atmospheric 
depth at which the shower attains its maximum size, X max . 
For a given shower, the position of X max will depend on 
the depth of the first interaction of the primary in the atmo- 
sphere and the depth that it takes the cascade to develop. 
Thereby, it will depend not only on the primary mass, but 
also on the cross section of the primary particle with air and 
on features of hadronic interactions at high energies. This 
important caveat should be kept in mind when discussing 
the mass composition of cosmic rays, i.e. interpretation of 
shower observables in terms of primary mass are subject 
to deficiencies of hadronic interaction models employed 
in air shower simulations. Besides the position of X max , 
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its fluctuations on shower-by-shower basis, RMS(X max ), 
show strong sensitivity to the primary mass. 

For the analysis, again hybrid data are used and the shower 
profiles are required to be good fits to a Gaisser-Hillas func- 
tion, as deviations could indicate the presence of residual 
clouds. Both (X max ) and its RMS show a characteris- 
tic change at E 5 • 10 18 eV indicating an increasingly 
heavier composition when compared to EAS simulations 
[3, P. Facal] (c.f. Fig. 3). It is well known that MC predic- 
tions are more uncertain for the (X maK ) than for the fluc- 
tuations. This is mainly due to the additional dependence 
of (X max ) on the multiplicity in hadronic interactions. In 
Fig. 2 we therefore compare the shape of the distributions, 
I max — (X max ), to MC predictions for different compo- 
sitions and hadronic interaction models. As can be seen, 
in this representation the various models predict a nearly 
universal shape. At low energy, the shape of the distribu- 
tion is compatible with a very light or mixed composition, 
whereas at high energies, the narrow shape would favor a 
significant fraction of nuclei (CNO or heavier). 

Fluorescence telescopes are the only observational tool cur- 
rently enabling direct measurements of the shower maxi- 
mum X max . Unfortunately, those data suffer from statistics 
because of the duty cycle being only ~ 15%. However, 
surface detectors, operated 24 hours a day, also provide ob- 
servables which are related to the longitudinal shower pro- 
file. These observables are subject to independent system- 
atic uncertainties (both experimentally and theoretically). 
The higher statistics allow us to extend these measurements 
to higher energies than possible with the FD. 

For each SD event, the water-Cherenkov detectors record 
their signals as a function of time. Since muons travel in 
almost straight lines whereas the electromagnetic particles 
suffer more multiple scattering on their way to ground, the 
first part of the signal is dominated by the muon compo- 
nent. Due to the absorption of the electromagnetic (EM) 
component, the number of these particles at the ground de- 
pends, for a given energy, on the distance to the shower 
maximum and therefore on the primary mass. In conse- 
quence, the time profile of particles reaching ground is sen- 
sitive to the cascade development: the higher the produc- 
tion height, the narrower the time pulse. The time distri- 
bution of the SD signal is characterised by means of the 
risetime, ti/%, which depends on the distance to the shower 
maximum, the zenith angle 6 and the distance to the core 
r. In a first step, the zenith angle at which the risetime 
asymmetries between the inner and outer SD stations of a 
shower become maximal is is analyzed. This angle O max 
is then, in a second step, related to the shower maximum 
using a subset of hybrid events [3, D. Garcia-Pinto]. Using 
this correlation it is possible to measure the shower evolu- 
tion with surface detector data, in a similar way to that in 
which the SD energy calibration is performed for a subset 
of events with the FD data. The result is shown in Fig. 3. 

Not only the risetime of signals in the SD tanks, but also 
the arrival time of particles with respect to the shower 
front plane contains information about the position of the 
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Figure 3: Results on shower evolution sensitive observables 
compared with models prediction. The error bars correspond to 
the statistical uncertainty. The systematic uncertainty is repre- 
sented by the shaded bands. 

shower maximum. A method for reconstructing the so- 
called Muon Production Depth (MPD), i.e. the depth at 
which a given muon is produced, measured parallel to the 
shower axis, using the FADC traces of detectors far from 
the core, has been presented in [4]. From the MPDs an 
observable can be defined, as the depth along the 

shower axis where the number of produced muons in a 
shower reaches a maximum. The method is currently re- 
stricted to inclined showers where muons dominate the sig- 
nal at ground level. Once the MPDs are obtained for each 
event, the value of X£ ax is found by fitting a Gaisser-Hillas 
function to the depth profile. The results of (X^ ax ) pre- 
sented in Fig. 3 are restricted to zenith angles between 55° 
and 65° and use timing information only for detectors far 
from the core (r > 1800 m). Because of this distance re- 
striction, the effective energy range for which the method 
can presently be applied is limited to E > 2 • 10 19 eV. 
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Figure 4: Proton-air cross section compared to other measure- 
ments and model predictions (for references see [3, R. Ulrich]). 
The inner error bars are statistical only, while the outer include all 
systematic uncertainties for a helium fraction of 25 % and 10 mb 
photon systematics. 

The measured values of (X£ ax ) are presented in the upper 
panel of Fig. 3. It is important to point out, that the predic- 
tions of X^ ax from different hadronic models would not 
be affected if a discrepancy between a model and data [3, 
J. Allen] is limited to the total number of muons. How- 
ever, differences in the muon energy and spatial distribu- 
tion would modify the predictions. 

With this caveat concerning hadronic interaction models, 
one might infer the primary composition from the data 
on the longitudinal air shower development presented in 
Fig. 3. The evolution of (X max ), 9 max and (X^ ax ) with 
energy is similar, despite the fact that the three analyses 
come from completely independent techniques that have 
different sources of systematic uncertainties. Concerning 
the RMS of X max , a variety of compositions can give rise 
to large values of the RMS, because the width of the X max 
is influenced by both, the shower-to-shower fluctuations of 
individual components and their relative displacement in 
terms of (X max ) [5]. However, within experimental un- 
certainties, the behaviour of (X max ), O max and (X^ ax ) as 
shown in Fig. 3 is compatible with the energy evolution of 
RMS(X max ). In particular, at the highest energies all four 
analyses show consistently that our data better resemble the 
simulations of heavier primaries than pure protons. 

4 p-air cross section and tests of hadronic in- 
teraction models 

One of the biggest challenges for a better understanding of 
the nature of ultra-high energy cosmic rays is to improve 
the modeling of hadronic interactions in air showers. None 
of the current models is able to consistently describe cos- 
mic ray data, which most importantly prevents a precise 
determination of the primary cosmic ray mass composition. 

Studies to exploit the sensitivity of cosmic ray data to the 
characteristics of hadronic interactions at energies beyond 
state-of-the-art accelerator technology began over 50 years 



ago. The property of interactions most directly linked 
to the development of extensive air showers is the cross- 
section for the production of hadronic particles (e.g. [6]). 
To reconstruct the proton-air cross-section based on hy- 
brid data, we analyse the shape of the distribution of the 
largest values of the depth of shower maxima, X max . This 
tail of the X max -distribution, that contains the 20 % of 
deepest showers, exhibits the expected exponential shape 
dN/dX max oc exp — (X max /A 2 o)- It is directly related to 
the p-air cross section via <7 p _ a i r = (m a i r )/A2o- In prac- 
tice, to properly account for shower fluctuations and detec- 
tor effects, the exponential tail is compared to Monte Carlo 
predictions. Any disagreement between data and predic- 
tions is then attributed to a modified value of the proton-air 
cross-section [3, R. Ulrich]. In this analysis, the energy in- 
terval is restricted to 10 18 to 10 18 5 eV which corresponds 
to a center-of-mass energy in the nucleon-nucleon system 
of y/s = 57 TeV. This interval has been chosen because 
of high statistics in the data and because of the composi- 
tion being compatible with a dominance of protons (see. 
Sec. 3). A possible contamination of He primaries could 
mimic a larger cross section (e.g. by 20 mb for 20% He 
contamination) while a photon contamination could reduce 
the cross section by at most lOmb. Combining the results 
one finds 

ffp-air = (505 ± 22 stat (tiDsyst) mb 

at a center-of-mass energy of 57 ± 6 TeV. This result is 
shown in comparison to other data and models in Fig. 4. 

The result favors a moderately slow rise of the cross 
section towards higher energies, well in line with recent 
results from LHC (e.g. [7]). A conversion of the derived 
Op-ah measurement into the more fundamental cross- 
section of proton-proton collisions using the Glauber 
framework [8, 9] will be published elsewhere. 

The importance of hadronic interaction models to measure- 
ments of the cosmic ray mass composition has been ad- 
dressed in Sect. 3. In particular, muons in extensive air 
showers are subject to large theoretical uncertainties due 
to our limited knowledge of multi-particle production in 
hadronic interactions. However, hybrid data can be used to 
constrain the models and to uncover deficiencies in describ- 
ing features of EAS data. When measuring the longitudinal 
profile (LP) in a golden hybrid event, we construct a li- 
brary of simulated air-shower events with the same shower 
geometry where the LP of each simulated event matches 
a measured one. The measured LP constrains the natural 
shower-to-shower fluctuations of the distribution of parti- 
cles at ground. This allows the ground signals of simulated 
events to be compared to the ground signals of measured 
events on an event-by-event basis. An example of such an 
analysis is shown in Fig. 5 [3, J. Allen]. Here, the LP of a 
measured event is compared to p and Fe simulations, each 
providing a good fit to the data. The bottom panel shows 
the corresponding signals in the SD. The ratio of the mea- 
sured signal at 1000m from the shower core, S'(IOOO), to 
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Figure 5: Top: A longitudinal profile measured for a hybrid event 
and matching simulations of two showers with proton and iron 
primaries. Bottom: A lateral distribution function determined for 
the same hybrid event as in the top panel and that of the two sim- 
ulated events. 

that predicted in simulations of showers with proton pri- 



maries. 



S(1000) sin 



is 1 .5 for vertical showers and grows to 
around 2 for inclined events. The ground signal of more- 
inclined events is muon-dominated. Therefore, the increase 
of the discrepancy with zenith angle suggests that there is 
a deficit of muons in the simulated showers compared to 
the data. The discrepancy exists for simulations of showers 
with iron primaries as well, which means that the ground 
signal cannot be explained only through composition. 

This finding is corroborated by direct estimations of the 
muon number from the signal traces in the SD as well 
as by making use of universality of the muonic to elec- 
tromagnetic signal S , A1 (1000)/5 em (1000) for fixed vertical 
depth of the shower [3, J. Allen]. Moreover, this purely- 
observational estimation of the muonic signal in data, is 
compatible with results obtained from inclined showers [3, 
G. Rodriguez]. This is best illustrated by Fig. 6. Here, A^g 
is defined as the ratio of the total number of muons, N^, 
in the shower with respect to the total number of muons 
St E — 10 EeV given by a 2-dim reference distribution, 
N 19 = N^E, 6)/N™ a *>(E = 10 EeV, 9) which accounts 
for the geomagnetic spatial deviation of muons at ground. 

Thus, all of the analyses show a significant deficit in the 
number of muons predicted by simulations with proton pri- 
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Figure 6: Fit of a calibration curve Nig = A(E/WEeV) B . The 
constants A and B are obtained using the maximum-likelihood 
method. The contours indicate the constant levels of the p.d.f. 
integrated over zenith angle, corresponding to 10, 50 and 90% 
of the maximum value [1, H. Dembinski]. Calibration curves for 
protons QGSJETII (dot line) and iron EPOS 1 .99 (dashed line) are 
shown for comparison. 

maries compared to data. This discrepancy cannot be ex- 
plained by the composition alone, although a heavy com- 
position could reduce the relative excess by up to 40%. 
The increased sophistication of the methods gives further 
weight to the conclusions that, at the current fluorescence 
energy scale, the number of muons in data is nearly twice 
that predicted by simulations of proton-induced showers. 
The possible zenith angle dependence of A^ el suggests that, 
in addition to the number, there may also be a discrepancy 
in the attenuation and lateral distribution of muons between 
the simulations and data. 

5 Update of photon and neutrino upper lim- 
its 

The search, and possibly study, of high energy photons and 
neutrinos is of interest for at least three reasons: (i) Top- 
Down models of UHECR origin [10] including topological 
defect or super-heavy dark matter models predict a signif- 
icant fraction of photons and neutrinos at the highest en- 
ergies, (ii) they would provide a smoking-gun signature of 
the GZK-effect because of the decay of charged and neu- 
tral pions created in photo-pion production, and (iii) they 
would open a new window to the most extreme Universe 
by possibly seeing point sources in the sky. A search for 
their signatures has thus been part of the Auger research 
program from the beginning. 

The search for EeV photons presented at this meeting is 
based on hybrid events. Due to the FD duty cycle the 
event statistics is reduced compared to the SD-only detec- 
tion mode. However, the hybrid detection technique pro- 
vides a precise geometry and energy determination with the 
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Figure 7: Upper limits on the photon flux above 1, 2, 3, 5 and 10 
EeV derived in this work (red arrows) compared to previous limits 
from Auger, from AGASA (A), and Yakutsk (Y). The shaded re- 
gion and the lines give the predictions for the GZK photon flux 
and for top-down models (TD, Z-Burst, SHDM and SHDM'). 
(See [13, M. Settimo] for references.) 

additional benefit of allowing to reduce the energy thresh- 
old for detection to about 10 18 eV. To improve the photon- 
hadron discrimination power over measurements of X max 
only, the differences in the lateral distribution functions for 
photons and hadrons measured by the SD have been con- 
sidered by analyzing the observable, 5&, defined in [1 1]. To 
reject misreconstructed profiles, only time periods with the 
sky not obscured by clouds, and with a reliable measure- 
ment of the vertical optical depth of aerosols, are selected. 
On the SD side we require at least 4 active stations within 
2 km of the hybrid reconstructed axis. This prevents an un- 
derestimation of Sb (which would mimic the behavior of a 
photon event) due to missing or temporarily inefficient de- 
tectors. For the classification of photon candidates, a Fisher 
analysis trained with a sample of a total of ^30000 pho- 
ton and proton CORSIKA [12] showers generated accord- 
ing to a power law spectrum between 10 17 and 10 20 eV is 
performed [13, M. Settimo]. The Fisher response distribu- 
tions for photon and proton primaries are well separated for 
all energies above 10 18 eV. Photon-like events in the data 
are then selected by applying an "a priori" cut to the upper 
50 % of the photon like events. This reduces the photon de- 
tection efficiency to 50 % but provides a conservative result 
in the upper limit calculation by reducing the dependence 
on the hadronic interaction models and on the mass com- 
position assumption. With this choice, the expected hadron 
contamination is about 1 % in the lowest energy interval 
(between 10 18 and 10 18 5 eV) and it becomes smaller for 
increasing energies [13, M. Settimo]. 

Applying the method to data, 6, 0, 0, and photon can- 
didates are found for energies above 1, 2, 3, 5 and 10 EeV. 
We checked with simulations that the observed number 
of photon candidates is consistent with the expectation 
for nuclear primaries, under the assumption of a mixed 
composition. The corresponding 95 % CL upper limits 
on the photon flux $9 5Ci integrated above an energy 
threshold Eq are shown in Fig. 7. To be conservative, a 



minimum value of the exposure above Eq is used and 
a possible nuclear background is not subtracted for the 
calculation of N® 5CL . The flux limits shown in Fig. 7 or 
likewise the derived limits on the photon fraction of 0.4%, 
0.5%, 1.0%, 2.6% and 8.9% above 1, 2, 3, 5 and 10 EeV, 
significantly improve previous results at the lower energies 
and rule out exotic models of UHECR origin, except for 
the Z-burst model of Ref. [14]. While the focus of the 
current analysis was the low EeV range, future work will 
be performed to improve the photon-hadron separation 
also at higher energies using further information provided 
by the SD. 

The surface detector is well suited also to search for ultra- 
high energy neutrinos in the sub-EeV energy range and 
above. Neutrinos of all flavours can interact in the atmo- 
sphere and induce inclined showers close to the ground 
(down-going). The sensitivity of the SD to tau neutrinos 
is further enhanced through the "Earth-skimming" mech- 
anism (up-going). Both types of neutrino interactions can 
be identified through the broad time structure of the signals 
induced in the SD stations. 

The analysis starts with the inclined shower selection 
(down-going: 9 > 75° and Earth-skimming 9 < 96°). 
These showers usually have elongated patterns on the 
ground along the azimuthal arrival direction. A length L 
and a width W are assigned to the pattern and a cut on 
their ratio L/W is applied. We also calculate the apparent 
speed V of an event using the times of signals at ground 
and the distances between stations projected onto L. Fi- 
nally, for down-going events, we reconstruct the zenith an- 
gle 9 ICC . After this pre-selection, the FADC signal-traces of 
the SD stations are analyzed to search for so-called "young 
showers" with a broad time structure. To optimize the dis- 
crimination power, again a Fisher discriminant method is 
used and trained to a subset of data. The identification ef- 
ficiency for the set of selection cuts applied to the data de- 
pends on the neutrino energy E v , the slant depth D from 
ground to the neutrino interaction point, the shower geom- 
etry, the neutrino flavour {v e , v^, or v T ), and is different for 
CC- and NC-type interactions, see [15]. 

Using the independent sets of identification criteria that 
were designed to search for down- and up-going neutri- 
nos in the data collected from 1 January 2004 to 31 May 
2010, no candidate was found [13, Y. Guardincerri]. As- 
suming a differential flux f(E v ) = kE~ 2 , we place a 
90% CL upper limit on the single flavour neutrino flux of 
k < 3.2 x 10 -8 GeV cm~ 2 s _1 sr _1 in the energy interval 
1.6 x 10 17 cV — 2.0 x 10 19 eV, based on Earth-skimming 
neutrinos and k < 1.7 x 10~ 7 GcV cm~ 2 s _1 sr _1 in the 
energy interval 1 x 10 17 eV — 1 x 10 20 eV, based on 
down-going neutrinos (see Fig. 8). The optimistic fluxes 
for p-primaries shown in this figure are accessible for the 
proposed lifetime of the Pierre Auger Observatory. The 
transition models and sources with a dominance of heavy 
primaries would be challenging to reach for any of the cur- 
rently operating experiments. 
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Figure 8: Differential upper limits (90% CL. per half decade of 
energy) from the Pierre Auger Observatory for a diffuse flux of 
down-going v (2 yr of full Auger) and Earth-skimming v T (3.5 
yr of full Auger [13, Y. Guardincerri]). Limits from other experi- 
ments and expected fluxes are also shown (see [5] for references). 

With no candidate events found in the search period, we can 
also place a limit on the UHE neutrino flux from a source at 
declination S. Since the sensitivity to UHEi^s is limited to 
large zenith angles, the rate of events from a point source in 
the sky depends strongly on its declination. In both Earth- 
skimming and down-going analyses the sensitivity yields a 
broad "plateau" spanning A6 ~ 110° in declination with 
the highest sensitivity reached at 5 ~ ±55°. The present 
flux limits do not yet allow us to constrain models of UHE^ 
production in the jets and the core of CenA [16]. 

6 Anisotropics 

One of the keys to understanding the nature of UHECRs is 
their distribution over the sky. This distribution depends 
on the location of the UHECR sources, as well on the 
UHECR mass composition and large-scale magnetic fields, 
both Galactic and extragalactic. Despite significant efforts, 
none of these issues is well understood at present. Obser- 
vation of the suppression of the CR flux at the highest ener- 
gies (c.f. Sect. 2) and its interpretation in terms of the GZK- 
effect suggests that the closest sources of UHECRs are sit- 
uated within the GZK volume of d G zK & 100 Mpc. At 
these scales the matter distribution in the Universe is inho- 
mogeneous, and so must be the distribution of the UHECR 
sources. If propagation of UHECRs at these distances is 
quasi-rectilinear, anisotropies would be expected, show- 
ing variations at large angular scales and possibly point 
sources. 

The Pierre Auger Collaboration in fact reported [17] di- 
rectional correlations of UHECR at E > 5.5 • 10 19 eV 
with AGN from the Veron-Cetty-Veron catalog [18] within 
75 Mpc on an angular scale of 3.1° at the 99 % CL. The 
optimal parameters were found using a exploratory scan 
(Period I) and independent data (Period II) showed 8 of 1 3 
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Figure 9: The most likely value of the degree of correlation 
Pdata = k/N is plotted as a function of the total number of 
time-ordered events (excluding those in period I). The 68%, 95% 
and 99.7% confidence level intervals around the most likely value 
are shaded. The horizontal dashed line shows the isotropic value 
Piso = 0.21 and the full line the current estimate of the signal 
Pdata = 0.33 ± 0.05. The black symbols show the correlation 
fractions bins of independent 10 consecutive events. 

events correlating. An update then yielded 21 of 55 events 
(Period II+III) correlating for the same parameter set. Here, 
we present the latest update including data up to June 201 1 
(c.f. Fig. 9) which yields a total of 28 of 84 events (Pe- 
riod II+III+IV) showing a correlation on a 3.1° -scale with 
a nearby AGN. The overall correlation strength thus de- 
creased from (69l 13 ) % initially to (33 ± 5)%. The chance 
probability of observing such a correlation from a random 
distribution remains below 1 %. Cumulative plots are of- 
ten misleading and Fig. 9 may be interpreted as a signal 
that is fading away. Thus, the superimposed black symbols 
show in addition the averages of 10 independent consec- 
utive events. Obviously, the first bin is an upwards fluc- 
tuation by about 3er from the mean of all events while the 
rest of the dataset does not show any peculiarity. Evidently, 
more data is needed to arrive at a definite conclusion. 

Interestingly, at this meeting the Telescope Array Collab- 
oration presented, though with a much lower exposure of 
~ 15% of Auger, an analysis of the northern sky adapting 
the parameters from the Auger collaboration (for a recent 
update see [19]). With 11 correlating events of 25 being 
above their energy threshold, they find a signal strength of 
44%. Correcting this for the larger chance probability of 
Piso = 0.24 compared to 0.21 in Auger, a good agreement 
of the data sets can be concluded. However, the TA events 
alone can originate from an isotropic distribution with a 
chance probability of about 2 %. The combined Auger and 
TA chance probability of observing such a correlation is at 
the 10" 3 level. 

The sky region around Centaurus A is populated by a larger 
number of high energy events compared to the rest of the 
sky, with the largest departure from isotropy at 24° around 
the center of Cen A with 19 events observed and 7.6 ex- 
pected for isotropy. However, a Kolmogorov-Smirnov test 
shows a chance probability for this to occur at a level of 
4%. Similarly, a search for directionally-aligned events 
(or 'multiplets') expected from sets of events coming from 
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the same source after having been deflected by interven- 
ing coherent magnetic fields shows one 12-plet with an en- 
ergy threshold of 20 Ee V. The probability that it appears by 
chance from an isotropic distribution of events is again 6 %. 
Thus, there is no significant evidence for the existence of 
correlated multiplets in the present data set [13, G. Golup] 
and [20]. It will be interesting to check if some of the ob- 
served multiplets grow significantly or if some new large 
multiplets appear. If one of them were a real multiplet, 
doubling the statistics should double its multiplicity, i.e. the 
significance does not increase as >/~N but much faster. 

The Pierre Auger Observatory has sensitivity also to neu- 
tron fluxes produced at cosmic ray acceleration sites in the 
Galaxy. Because of relativistic time dilation, the neutron 
mean decay length is (9.2 x E) kpc, where E is the neu- 
tron energy in EeV. A blind search over the field of view 
of the Auger Observatory for a point-like excess yields no 
statistically significant candidates. The galactic center is a 
particularly interesting target because of the presence of a 
massive black hole. The results for the window centered on 
it and for E > 1 EeV shows no excess with a 95% CL up- 
per limit on the flux from a point source in this direction of 
0.01 km -2 yr" 1 [13, B. Rouille d'Orfeuil], which updates 
the bounds obtained previously [21]. We note that for di- 
rections along the Galactic plane the upper limits are below 
0.024 km~ 2 yr~\ 0.014 km~ 2 yr" 1 and 0.026 km~ 2 yr" 1 
for the energy bins [1 - 2] EeV, [2 - 3] EeV and E > 1 EeV, 
respectively. 

A targeted search has also been performed to test potential 
sources of galactic cosmic rays, such as SNR, pulsars and 
Pulsar Wind Nebulae (PWN). The candidate sources are 
expected to be strong gamma-ray emitters at GeV and TeV 
energies. For this reason, we apply a neutron search also 
to Galactic gamma-ray sources extracted from the Fermi 
LAT Point Source Catalog [22] and the H.E.S.S. Source 
Catalog 1 . Targets were selected among the sources located 
in the portion of the Galactic plane, defined as | fo| < 10°, 
covered by the FOV of the SD, and located at a distance 
shorter than 9 kpc (A n at 1 EeV) [13, B. Rouille d'Orfeuil]. 
The stacked signal computed from the SD data at the 
positions of the two sets of sources under study and for the 
same energy bins used in the Galactic plane search has not 
yet yielded an excess. 

Besides searching for point sources of charged cosmic rays 
or neutrons, the large scale distribution of arrival directions 
of CRs represents another important tool for understanding 
their origin. Using data from the SD array, upper limits 
below 2% at 99% CL. have been recently reported for 
EeV energies on the dipole component in the equatorial 
plane [23]. Such upper limits are sensible, because cos- 
mic rays of galactic origin, while escaping from the galaxy 
in this energy range, might generate a dipolar large-scale 
anisotropy with an amplitude at the % level as seen from 
the Earth. Even for isotropic extragalactic cosmic rays, a 
large scale anisotropy may remain due to the motion of our 
galaxy with respect to the frame of extragalactic isotropy. 
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Figure 10: Upper limits on the anisotropy : equatorial dipole 
component d± as a function of energy from Auger. Results from 
EAS-TOP, AGASA, KASCADE and KASCADE-Grande exper- 
iments are also displayed, in addition to several predictions (see 
[13, H. Lyberis] for references. 

This anisotropy would be dipolar in a similar way to the 
Compton-Getting effect [24] in the absence of the galactic 
magnetic field. 

An update of the results of searches for first harmonic mod- 
ulations in the right ascension distribution of cosmic rays 
is presented in Fig. 10 [13, H. Lyberis]. The upper lim- 
its at 99% CL obtained here provide the most stringent 
bounds at present above 2.5 x 10 17 eV. Some predictions 
for anisotropies arising from models of both galactic and 
extragalactic cosmic ray origin are included in the plot to- 
gether with data from other experiments. In models A and 
S (A and S standing for 2 different galactic magnetic field 
symmetries [25]), the anisotropy is caused by drift motions 
due to the regular component of the galactic magnetic field, 
while in model Gal [26], the anisotropy is caused by purely 
diffusive motions due to the turbulent component of the 
field. Some of these amplitudes are challenged by our cur- 
rent sensitivity. For extragalactic cosmic rays considered 
in model C-G Xgal [27], the motion of our galaxy with re- 
spect to the CMB (supposed to be the frame of extragalactic 
isotropy) induces the small dipolar anisotropy (neglecting 
the effect of the galactic magnetic field). 

While the measurements of the amplitudes do not provide 
any evidence for anisotropy, it is interesting to note that 
the phase shown in Fig. 1 1 suggests a smooth transition be- 
tween a common phase of ~ 270° below 1 EeV and an- 
other phase (right ascension ~ 100°) above 5 EeV. This 
is potentially interesting, because, with a real underlying 
anisotropy, a consistency of the phase measurements in or- 
dered energy intervals is indeed expected with lower statis- 
tics than that required for the amplitudes to significantly 
stand out of the background noise. Applying a Likelihood 
test leads to a probability of ~ 10 -3 of observing this from 
a random distribution. However, since we did not perform 
an a priori search for such a smooth transition in the phase 

1 . http://www.mpi-hd.mpg.de/hfm/HESS/pages/home/sources/ 
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Figure 1 1 : Phase of the first harmonic as a function of energy. 
The dashed line, resulting from an empirical fit, is used in the 
likelihood ratio test (see text) [13, H. Lyberis], 



measurements, no confidence level can be derived from this 
result. 

The infill surface detector array which is now operating at 
the Pierre Auger Observatory will allow us to extend this 
search for large scale anisotropics to lower energy thresh- 
olds. 
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Figure 12: Ten seconds average of the Auger scaler rate for the 
27 Feb 2010 Chile major 8.8 magnitude earthquake. A strong 
24 a decrease is found 90 ± 2 (stat) seconds afterwards, compat- 
ible with the time delay expected for seismic S-waves traversing 
the distance from the epicentre to the Auger Observatory. 
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7 Serendipity observations and interdisci- 
plinary science 

The hybrid character of the Auger Observatory, but also 
the surface and fluorescence detectors themselves allow a 
number of studies beyond cosmic ray physics. 

A first study, using the so-called scaler mode data of the 
SD has been presented in [28] and is updated at this con- 
ference [13, H. Asorey]. It uses the count rates of low 
energy secondary cosmic ray particles (deposited energy 
> 15MeV and 15 < E dcp < 100 Me V, using two differ- 
ent triggers) recorded continuously for self-calibration pur- 
poses for all of the 1660 SD stations. With each detector 
recording about 3600 Hz, we thus record a total of about 
6 MHz so that even very small changes of the rates due to 
atmospheric and solar changes can be monitored. This en- 
ables the SD of Auger to address questions of solar cosmic 
rays and allows to study Forbush events. A good agree- 
ment between neutron monitor and the scaler data is found 
when accounting for different geomagnetic cut-offs of de- 
tectors located at different latitudes and for different effec- 
tive energy thresholds of neutron monitors and the Auger 
SD stations [28]. 

Instead of using averaged scaler rates for the whole array, 
it is also possible to study the scaler rate of individual sta- 
tions, in order to study the propagation of some phenomena 
across the Auger SD, like the crossing of a storm over the 
3000 km 2 of the array. This is because the flux of secondary 
particles changes as the pressure front moves from across 
the detector field. Additional analyses to study the influ- 
ence of the variation of electric fields on the flux of EAS 
particles are currently being carried out as well. 
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Figure 13: FD camera image for 4 consecutive time windows as 
indicated. It shows the time evolution of an elve located at about 
80 km altitude at a distance of 580 km from the observatory [29, 
A. Tonachini]. 

Interestingly, the 8.8 magnitude earthquake in Chile on 
27 Feb 2010 06h34 UTC with the epicentre located about 
300 km SW from the Auger Observatory left traces in the 
SD as well. The averaged scaler rate for the whole ar- 
ray and also for individual stations showed a 24 a decrease 
beginning (90 ± 2) seconds after the earthquake. This de- 
lay is compatible with the propagation of seismic S-waves 
over that distance. The scaler rate from 6hl5 to 6h45 UTC 
is shown in figure 12. Although other minor quakes have 
been recorded by seismographs near the SD, no other sim- 
ilar effects have been found in 6 years of data. Detailed 
analyses to identify the causes of the observed drop in 
the scaler rate are underway. These include simulations 
and shaking tests of selected detectors in the array. Af- 
ter 6 hours, the scaler rate recovered to the mean value for 
February 2010. 
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Also quite unexpectedly, during a normal FD data taking 
shift an unusual event has been observed with a well de- 
fined space-time structure: a luminous ring starting from 
a cluster of pixels, and expanding in all directions [29, 
A. Tonachini]. Usually, such kind of events lasting for 
much longer than 70 /is and with such a high multiplicity 
are rejected by the T2 trigger because of being caused by 
lightning with high probability. Due to this rejection, only 
three of such unusual events were recorded. By careful re- 
construction of the timing, these events could be identified 
as elves originating from lightning in the western part of 
Argentinia. Elves are transient luminous phenomena origi- 
nating in the D-layer of the ionosphere, high above thunder- 
storm clouds, at an altitude of approximately 90 km. With 
a time resolution of 100 ns and a space resolution of about 
1 degree, the FD can provide an accurate 3D measurement 
of elves for thunderstorms which are below the horizon. To 
improve the detection efficiency for such kind of interest- 
ing and not well understood phenomena, a dedicated trig- 
ger will be implemented in the future. 

The Auger Observatory allows us to perform a number of 
further interdisciplinary science studies, mostly related to 
atmospheric sciences (study of aerosols, atmospheric grav- 
ity waves, etc.) but includes also biological studies in the 
pampa as well as related studies of earthquakes either di- 
rectly by its instrumentation operated or indirectly by pro- 
viding infrastructure for non-cosmic ray scientific commu- 
nities. 

8 Summary and Conclusions 

The Pierre Auger Observatory has reached a cumulative ex- 
posure of more than 25 000 km 2 sr yr by the time of writing 
this article. This exceeds by far the total statistics recorded 
by all other observatories. A great deal of new insights are 
provided by these data, but many new questions have ap- 
peared. This is most prominently about the origin of the 
suppression of the CR flux at highest energies and - related 
to this - the mass composition and anisotropies at the high- 
est energies. The Observatory will continue to collect data 
with unprecedented precision for several more years and 
it is hoped that these data will help to unravel the puzzles 
about the most energetic particles in nature. 
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